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Abstract: We calculate the total cross section for deep inelastic scattering (DIS) on a nucleus at 
high energy for a strongly coupled M = 4 super Yang-Mills theory using AdS/CFT correspondence. 
r^, ; In analogy to the small coupling case we argue that at high energy the total DIS cross section is 
^ ; related to the expectation value of the Wilson loop formed by the quark-antiquark dipole. We 
model the nucleus by a metric of a shock wave in AdS 5 . We then calculate the expectation value of 
the Wilson loop by finding the extrema of the Nambu-Goto action for an open string attached to 
the quark and antiquark lines of the loop in the background of an AdSs shock wave. We find three 
extrema of the Nambu-Goto action: the string coordinates at the extrema are complex-valued and 
are given by three different branches of the solution of a cubic equation. The physically meaningful 
solutions for the total DIS cross section are given either by the only branch with a purely imaginary 
string coordinate in the bulk or by a superposition of the two other branches. For both solutions 
we obtain the forward scattering amplitude iV for the quark dipole-nucleus scattering. We study 
the onset of unitarity with increasing center-of-mass energy and transverse size of the dipole: we 
observe that for both solutions the saturation scale, while energy-dependent at lower energies, at 
very high energy becomes independent of energy/Bjorken-x. The saturation scale depends very 
strongly on the atomic number of the nucleus as Q s ~ A 1 / 3 . 
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1. Introduction 



Over the past years there has been some significant progress made in our understanding of deep 
inelastic scattering (DIS) at high energy/small Bjorken-x in the framework of the Color Glass Con- 
densate (CGC) @, 1, & § 1, g, g 1 1, H, |n], 0J, ^ 0, H, 0, 0, H H H P, H, H, H. It 
was shown that at moderately small XBj the total DIS cross section is well- described in the quasi- 
classical (Glauber-Mueller (GM)/McLerran-Venugopalan (MV)) [^, [|, [|, U approximation includ- 
ing all multiple rescatterings each taken at the lowest (two-gluon) order. At smaller xsj the energy- 
dependence of the cross section comes in through the non- linear small-x^j Balitsky-Kovchegov (BK) 
20l [21], [lj], [18| and Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) 



10| , PI , |r2[ [13| , [14] , |15| , pT6|| quantum evolution equations. The BK and JIMWLK evolution equa- 

evolution equation at 



tions unitarize the linear Balitsky-Fadin-Kuraev-Lipatov (BFKL) [[25 
high energies in the large- N c limit (BK) and beyond (JIMWLK). 

The main principle of the CGC approach is the existence of a momentum scale Q s , called the 
saturation scale, which regulates the infrared divergences in the total DIS cross section and in other 

The saturation scale 



relevant observables for high energy collisions (for a review see [[22], [24] 
grows as a power of the center-of-mass energy of the collision (a power of Bjorken x) as follows from 
BK and JIMWLK equations, and as a power of the atomic number A for DIS on the nucleus, as can 
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be seen in the quasi-classical GM and MV approaches. More specifically, for small- x Bj evolution at 
the leading logarithmic (a s lnl/x^j) approximation (LLA) at fixed coupling Q 2 ~ (l/xBj) Sas A 1 / 3 , 
where a s is the strong coupling constant and 5 ~ 4.8 |28|, |29|| . Thus at high energies and/or for large 
nuclei the saturation scale can be large, much larger than Aqcd, cutting off the non-perturbative 
infrared (IR) effects. This justifies the perturbative approach and makes CGC physics perturbative. 

While the picture described above is theoretically robust and self-consistent, the physics of 
CGC has not yet reached the level of precision required to make unambiguous phenomeno logical 
predictions. The fixed coupling LLA approaches, while giving interesting qualitative predictions 
confirmed by HERA and RHIC data 



36|| , require many additional assumptions 



30|,^1|, |32|,g3j£4j, 

to describe the data quantitatively. At the same time it is known that NLO corrections to the BFKL 
kernel (and therefore to BK and JIMWLK kernels as well) are large [^, [38], 39| . They come in 



with a sign opposite to that of the LLA kernel, leading to cross sections decreasing with energy 
and other undesirable consequences [^0, [II] . It is believed that resummation of all higher order 
corrections should remedy this problem. This idea is supported by the success of resummation of 
collinear singularities in the all-order BFKL kernel performed in 



Recent calculation of the 



running coupling corrections to the BFKL, BK and JIMWLK |43|, [44], [|5|, (S| and the resulting 
successful phenomenology both for DIS and for heavy ion collisions at RHIC [|3(| f|7| also support 
the possibility that resummation of all higher order corrections to BK and JIMWLK evolution 
equations would improve the agreement with the data and resolve theoretical problems posed by 
NLO BFKL calculation. 

While it is not clear at all how to resum higher order corrections to the BK and JIMWLK 
kernels to all orders in QCD, one turns for guidance to other QCD-like theories, such as Af = 4 
super Yang-Mills (SYM), where one can perform calculations in the non-perturbative limit of large 
't Hooft coupling due to the Anti-de Sitter space/conformal field theory (AdS/CFT) correspondence 
p8| , f49[ , |50| , [5l| . The original work in this direction was performed by Janik and Peschanski ]52 
(see also pi M, EE 



)6[), who showed that at large 't Hooft coupling the scattering cross section 
mediated by a single pomeron exchange, corresponding to a single graviton exchange in the bulk 
of AdSs, grows as the first power of energy, cr to t ~ s. Comparing this result to the LLA BFKL 
prediction at small coupling of a to t ~ s ^ ln2 [p6 , 27], which is the same for both QCD and Af = 4 



SYM being entirely due to gluon dynamics, one concludes that the power of energy changes from 
to 1 as 't Hooft coupling A = An a s N c goes from being small to being large. The result of j52| 
was further expanded in [f57| . It also agrees with the extrapolation of the results of weak coupling 
resummations ||58j| . The conclusion one may draw from these results is that the cross sections should 
grow fastest with energy when the coupling is large. However, as in QCD the coupling is large at 
low momentum transfer Q 2 , one would conclude that the cross section should grow steeper with s 
as Q 2 gets lower. Such conclusion would completely contradict the existing DIS and proton-proton 
collisions phenomenology ||59|| , which shows that the exact opposite is true: the growth of cross 



sections with energy slows down with decreasing Q 2 . 
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We believe the resolution of this puzzle may lie in the physics of parton saturation/CGC. As 
one goes towards lower Q 2 , single-pomeron exchange approximation becomes invalid, and one has to 
resum multiple pomeron exchanges. These exchanges unitarize the total cross section, significantly 
reducing its growth with energy [F24]] . Hence to make the strongly-coupled dynamics consistent with 
phenomenology, saturation/CGC effects need to be included in the picture as well. The first steps 
in this direction were performed in pD| , |5T| , where the authors studied DIS on a static M = 4 



SYM at strong coupling by calculating a correlator of two .R-currents. Here we will present a 
different approach to DIS in AdS. We will argue below that the total DIS cross section at high 
energy can be related to an expectation value of the Wilson loop formed by the propagators of a 
quark and an anti-quark arising in the splitting of the incoming virtual photon. The expectation 
value of a Wilson loop can be found at strong coupling using the methods outlined in [|2|] (see also 



|, ^4], 0). We will model the nucleus by a shock wave metric in AdSs. The expectation value of 
the Wilson loop in the shock wave background is related to value of the Nambu-Goto action at the 
extremal string world-sheet [B2| for an open string connecting the quark and the anti-quark lines. 
We therefore extremize the open string action. Our calculation is similar to those performed in 
f66] , |67| , |68| , m [70], [71], |72|, [73], [74], [7^1 . We obtain the expectation value of the Wilson loop, which 
is easily related to the forward scattering amplitude N(r, Y) of a qq dipole of transverse size r on 
the nuclear target with the total rapidity interval of Y ~ In s ~ In 1/xBj- The obtained N(r, Y) is 
unitary and exhibits saturation transition when the dipole transverse size becomes sufficiently large. 
We extract the saturation scale from the obtained expression and notice that at very high scattering 
energy s, corresponding to large Y or small Bjorken-x, the saturation scale becomes independent 
of energy, contrary to the perturbative behavior outlined above. 

The paper is structured as follows. In Sect. ^ we discuss the relation of the total DIS cross 
section to the expectation value of a Wilson loop. We set up the problem in AdS5 space and 
introduce the shock wave metric. We look for the extrema of the open string connecting the quark 
and the anti-quark in the shock wave background in Sect. || We argue in Sect. |3| that for broad 
enough shock waves the Nambu-Goto action is dominated by static string configurations. Such 
approximation is justified for scattering on nuclei with sufficiently large atomic number A. We 
solve the static equations of motion for the open string and find that there are three possible 
complex-valued solutions for string coordinates extremizing Nambu-Goto action in the presence of 
the shock wave. Similar results are usually obtained in the quasi-classical approximation in quantum 
mechanics (see Ch. 131 in []7Efl). In the context of AdS/CFT similar properties have been found in 

We study the forward scattering amplitudes N(r, Y) (or, equivalently, the S-matrices S(r, Y)), 
resulting from these solutions in Sect. We notice that only one branch of the solution gives 
a physically meaningful unitary non-negative N(r, Y) for all r and Y. The branch has purely 
imaginary string coordinates in the bulk. The corresponding N(r, Y) is given by Eq. (|4.14|) and is 
plotted in Fig. [7[ We notice that a slight modification of Maldacena's prescription for the calculation 
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of Wilson loops shown in Eq. ( 4.29|) allows for a combination of two other branches to also give a 
meaningful and physical amplitude N(r,Y). The resulting amplitude in shown in Eq. ( |4.30| ) and 
is plotted in Fig. [|. 

The saturation scale Q s as a function of center-of-mass energy s is found in Sect. ^ as well and 
is shown in Figs. § and [10] for both solutions. For both solutions it appears that at intermediate 
energies the saturation scale grows with s. However the growth of Q s with s stops at very high 
energies, reminiscent of "saturation of saturation" conjectured in ||78|1 . We try to interpret this 
result in Sect. ^, where we also plot the resulting behavior of the saturation scale in the strong 
coupling sector of QCD in Fig. [II]. 



2. Setting up the Problem 



In DIS at high energies, when viewed in the target rest frame, the incoming virtual photon splits 
into a quark-antiquark pair, which then scatters on the proton or nuclear target |79|, |25|, |8(| . The 
process is depicted in Fig. [3]. The scattering of a qq pair on the target is eikonal at high energies 
||79| . p5| , 80, nj. The light-cone lifetime of the qq pair is much longer than the nuclear radius. Using 
light-cone perturbation theory |8l]. 82 one can therefore decompose the total scattering cross section 
of DIS into a convolution of a light-cone wave function squared ($) for a virtual photon to decay 
into a qq pair and the imaginary part of the forward scattering amplitude for the qq pair-target 
interaction (N) 0, g|, : 



d 2 r da 
2tt 



$(r,a,Q 2 ) d 2 b N(r,b,Y). 



(2.1) 



Here N(r, b, Y) is the imaginary part of the forward scattering amplitude for the scattering of a 
quark-antiquark dipole of transverse size r at center-of-mass impact parameter b on a target, where 
the total rapidity of the scattering process is Y = In 1/xsj with xbj the Bjorken x variable. a) 
is the light-cone wave function squared of the virtual photon with virtuality Q 2 splitting into a qq 
pair of transverse size r with the quark carrying a fraction a of the virtual photon's longitudinal 
(plus) momentum. (Boldface notation denotes two-dimensional vectors in the transverse plane.) 

It is important to stress that the factorization of Eq. fl2.1|) is independent of the strong interac- 
tion dynamics. The wave function squared $(r, a, Q 2 ) is purely due to electromagnetic interactions, 
and is very well known J75], |2"4]| . All the strong interaction dynamics is contained in N(r, b, Y). 

In perturbation theory, at the leading logarithmic level in a s lnl/x^, the scattering between 
the qq pair and the target is eikonal. One can therefore relate N(r, b, Y) to the expectation value 
of a fundamental Wilson loop |]| ^ ^ U, 



N{r,b,Y) = 1 - S(r,b,Y) 



(2.2) 
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photon 



Figure 1: The space-time picture of the forward amplitude for the deep inelastic scattering (DIS) at high 
energy. The collision axis is labeled x 5 , the time is x°. The virtual photon (7*) travels along one light 
cone and splits into a Wilson loop. The qq pair in the loop is separated only in the transverse plane: the 
separation shown along the x 3 -direction is for illustration purposes only. The Wilson loop then scatters 
on the target, traveling along the other light cone. As it is a forward amplitude, qq pair recombines back 
into 7* after the scattering. 



with the (real part of the) S-matrix 



S{r, b,Y) = ^Be(w(b+±r,b- \r, y\ 



(2.3) 



W (b + 2?*, b — gt*, Y) denotes a Wilson loop formed out of a quark line at b+ \r and an antiquark 
line at b — |r with the links connecting the two lines at plus and minus infinities on the light cone. 
In the frame where the target is at rest, the Wilson loop is oriented close to the light cone, with 
rapidity Y reflecting the "angle" with respect to the light cone p"9L P21j . The averaging in Eq. (127 



denoted by (...), is performed over all possible wave functions of the target. 

Neglecting the gauge links at the light-cone infinities by choosing an appropriate gauge, in QCD 
perturbation theory the Wilson loop can be rewritten in terms of Wilson lines [§, 10, 11, 12, [H| |14| , 
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m rra, ra, pi 



W (x, y, Y) = tr[U(x, Y) U\y, Y)\ (2.4) 



U(x,Y) = Pexpl -ig / dx^A" \ . (2.5) 



where the Wilson line is 



The integration contour runs from — oo to +00 close to the light cone with the "angle" defined by 
rapidity Y and being the gluon field of the target. The trace in Eq. ( |2.4| ) is in the fundamental 
representation. In perturbative LLA, in the absence of a target one has A 1 = 0, leading to U — 1, 
W = N c giving S = 1 and N = 0, indicating zero interaction cross section as expected. In the 
opposite limit of very strong interactions one should get the black disk limit, which corresponds to 
N — 1 and 5 = Q. Hence < N, S < 1. 

The QCD scattering problem as formulated by Eqs. (|2.2|) , (|273|) and ( |2.1|) and shown in Fig. [I] 
can be easily generalized to M = 4 SYM theories at strong coupling by means of the AdS / CFT 



prescription outlined in |62], Q (for a review see fl6"5]). Following [|62]| , we represent the quark 



and the anti-quark in the J\f=4 SYM theory as W-bosons coming from the breaking U(N C + 1) — > 
U (N c ) xU(l). Our "quarks" are going to be very massive and would not recoil when interacting with 
the nucleus, thus justifying the Wilson loop approximation for the forward scattering amplitude.* 
The gravity dual description of the A/" =4 SYM qq dipole corresponds to an open string in the AdSs 
space whose endpoints, the quark and the anti-quark, are located at the boundary of the AdSs 
space. Parameterizing the two-dimensional world sheet of the string by the coordinates (r, a), the 
location of the string in the five-dimensional world is given by 

X^ = X^(r,a), /i = 0,...,4. (2.6) 

The string's Nambu-Goto action is 



Sng = 7; 7 / drd(j^-detG af3 (2.7) 

2 71 a J 

with 

G aP = g^(X) d a X» dp X", a, (3 = a, r, (2.8) 



*Whilc the infinite mass of the "quarks" justifies the recoilless approximation for their interaction with the target, 
the question arises of applicability of this model to the description of real DIS processes with light quarks, which 
indeed recoil during the interaction. For instance, while at small coupling in the leading logarithmic a s lnl/xs^ 
approximation the recoilless approximation is justified pS, 21, Ya its validity becomes less clear at the subleading 



logarithmic order (order a 2 s In 1/xgj) |39|. We note that the recoil of the quarks at high energy only affects the impact 
factors (see e.g. j84|), and not the part of the interaction described by the small- x Bj evolution. Thus the calculation 
below, when applied to light quarks, should be understood as calculation of the evolution without the impact factors. 
While the impact factors are likely to be numerically important, we will proceed under the assumption that they 
will not affect the qualitative features of the obtained scattering amplitude. 
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where g^ v will be taken below to be the metric of the AdSs space in the presence of the shock wave 
and a' is the slope parameter. The AdS/CFT correspondence prescription of |6^] dictates that in 
the limit of large 't Hooft coupling 

\ = g\ M N c (2.9) 

with gyu the Yang-Mills coupling constant, and in the large- N c limit, the expectation value of the 
Wilson loop is given by the classical Nambu-Goto action of the open string in AdSs space 

(W c ) ~ e iSNG , (2.10) 

where C represents the closed contour spanned by the quark and the anti-quark at the string 
endpoints. 




target 



Figure 2: The same DIS process as in Fig. ffl, only viewed in the rest frame of the dipole, as considered 
in the text. Separation of the quark and the anti-quark lines in x 3 direction was put in the figure only to 
guide the eye: in the text both the quark and the anti-quark are located at the same x 3 . 

To complete the gravity dual description of the DIS scattering process we still have to determine 
the metric of the AdSs space in the presence of the nucleus. To that end we first introduce light-cone 
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coordinates 



x 



± 



X° ± X 3 

~7T 



(2.11) 



with x° the time direction. Following Janik and Peschanski ||85|| , we consider the ultra-relativistic 
nucleus as a shock wave moving along the positive x 3 direction with its energy momentum tensor 
given by 



(T__> 



N 2 

), 



(2.12) 



with all the other components being equal to zero. The physical interpretation of the scale /z, that 
has dimensions of mass cubed, and of the factor N 2 in Eq. ( |2.12|) will be discussed below. Then, 
according to holographic renormalization |86 |, and using Fefferman-Graham coordinates |[7|, the 
gravity dual description of the CFT energy-momentum tensor in Eq. (|2.12|) is given by the following 
line element 



ds 2 



L'- 



\ [ —2dx + dx + fj,5(x ) z 4 dx 2 + dx 2 ± + dz 2 ~\ 



(2.13) 



where dx\ = (dx 1 ) 2 + (dx 2 ) 2 with x l and x 2 being the coordinates transverse to the direction of 
motion of the nucleus, z is the coordinate describing the 5th dimension and L is the curvature 
radius of the AdS 5 space. In these coordinates the boundary of the AdS 5 space, and therefore the 
endpoints of the string representing the dipole, is at z = 0. 

From now on we will work in the qq dipole's rest frame. There the Wilson loop in Eq. ( |2.3j ) 
becomes time-like, as shown in Fig. |2|. We choose the dipole orientation to be parallel to the plane 
of the shock wave. We believe such configuration is the most relevant for studying DIS. Our choice 
of coordinates is such that the quark and anti-quark are located at x 2 = x 3 = z = and x\ = ± r/2 
respectively, with r the dipole's transverse size. We recall that, due to the eikonal approximation, 
the quark and the anti-quark remain fixed at those spatial locations throughout the collision process, 
thereby providing the appropriate Dirichlet boundary conditions for the string motion. Therefore, 
we find it convenient to parametrize the string coordinates ( |2.6| ) in terms of the proper time in the 
dipole rest frame, t, and x 1 = x in the following way: 



t 



a 



x G 



r r 
*2' 2J 



(2.14) 



a) -> (X° = t, X 1 = x, X 2 = 0, X 3 = 0, X 4 = z(t, x)) . (2.15) 

We note that the symmetries of the problem allow us to set X 2 = for all times, since nothing 
in the problem depends on x 2 . Analogous symmetry arguments indicate that the string coordinate 
along the collision axis, X 3 , is equal to zero before the collision with the target. A priori, there is no 
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dynamical argument allowing us to set X 3 = at all times. Indeed, the shock wave metric ( 2.1 3| ) is 
strongly singular at t = x 3 and, therefore, it is conceivable that after being hit by the nucleus the 
string develops some kind of motion in x 3 -direction, such as recoil or oscillation. However, it can 
be proven that for the particular dynamical limit that we are interested here, namely in the static 
limit that we shall discuss in detail in Section |3|, the dynamics of the X3 component of the string 
coordinates is trivial, allowing us to set it to a constant value, which we choose here to be equal to 
zero. In the general time-dependent case we would not expect X 3 to be equal to zero. The collision 
process as described in the previous paragraphs is represented in Fig. ||. 

Before proceeding any further let us pause to interpret our choices Eq. (|2.12|) and Eq. (|2.13| ), 
which provide the embedding of the ultra-relativistic nucleus into the AdS 5 space. First of all, as it 
was shown in |85 |, the shock wave metric in Eq. ( |2.13| ) is an exact solution of the Einstein equations 
in five dimensions with negative cosmological constant A c = — 6/L 2 , and, therefore, provides a 
plausible gravity representation of an ultra-relativistic nucleus. Secondly, let us relate the parameter 
fi in Eq. ( |2.12| ) to some physical quantities. The energy-momentum tensor in Eq. ( |2.12| ) corresponds 
to a nucleus made out of many massless ultra-relativistic point particles. Thus, we envision the 
nucleus as a set of A nucleons, each of them consisting of N 2 valence point-like charges moving with 
light-cone momentum p + along the positive x + direction. In our model the nucleus is infinite and 
homogeneous in the (x 1 , x 2 )-plane transverse to its direction of motion (x 3 ). The energy-momentum 
tensor of an ultra-relativistic point particle moving with momentum p + along the light cone at 

transverse position y and at the other light cone coordinate x~ = is T (x) = p + 5 2 (x — y) 5(x~). 

Summing over A" 2 gluons in each nucleon and over all A nucleons in the nucleus and averaging over 
transverse coordinates y we write 

A N 2 f AN 2 N 2 

{T nucleu S) = / d 2 yp+6 2 {x _ y)5{x - ) = = -JLA 1 ** A 2 P + 5(x~) (2.16) 

where S± ~ A 2 ! 3 tc R 2 n is the transverse area of the nucleus, Rn is the nucleon radius and A 2 = 

.2 /_ D 2 



2~k /ttR n is a characteristic transverse momentum scale. Comparing Eqs. ( |2.12| ) and ( |2.16| ) one 



can now identify the scale \i in Eq. fl2.12p in terms of the physical scales characterizing the nucleus 

ji = p + A 2 A 113 . (2.17) 

It is worth noting that if we had not assumed the existence of A" 2 valence particles per nucleon, 
then the perturbation to the metric of the empty AdSs space due to the presence of the nucleus 
would be suppressed by a factor of l/N 2 , becoming negligible in the strict N c — > oo limit. 

Eq. Q2.17D will serve as our expression for fi. A more general treatment for the energy-momentum 
tensor of a nucleus in the context of high-energy scattering has been recently proposed by the 
authors in [Q. In particular one could argue that the energy-momentum tensor of the nucleus at 
large coupling should not be limited to only the valence quarks (gluons) contribution. The strong 
gluon fields should also contribute to the energy- momentum tensor, and their contribution is likely 
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to be dominant at large coupling. The energy-momentum tensor of such fields is fixed by conformal 
invariance of Af = 4 SYM theory at hand. The structure of the energy-momentum tensor is the 
same as in an Abelian theory: it was calculated in Appendix A of ||88|| . Here we will simply argue 



that such energy-momentum tensor could be approximated by Eq. ( 2.12 ) with \i given by Eq. ( p. 17 ) 
times a factor of y/X. Such tensor would not qualitatively change our discussion below, but it will 
modify the expression ( |2.17] ) for u. that we will use in what follows. 

Under the approximation of a homogeneous nucleus having an infinite transverse extent the 
nuclear energy-momentum tensor is independent of the transverse coordinates, as reflected in 
Eq. ( |2.12|) . This means that all the relevant dynamical quantities become rotationally invariant and 
independent on the impact parameter of the collision, depending only on the dipole size, r = \r\, 
and on the collision energy: 



W (x, y,Y) -> W (r, Y) , N(r, b, Y) -> N(r, Y) , S{r, b, Y) -> S(r, Y). 



(2.18) 



Let us generalize the perturbative formulation of the scattering problem in Eq. 
large coupling case. Using Eqs. (|2.2j), ( |2-3p and ( |2.1(J| ) we write 



S(r,Y) = Re 



and 



to the 



(2.19) 



N{r,Y) = l-S{r,Y). 



(2.20) 



In Eq. ( |2.19D (W)^ is the static Wilson loop pictured in Fig. ^] and Sng(^) is the corresponding 
classical Nambu-Goto action of the open string evaluated in the background of the shock wave. 
(W)^q and Sng(h ~~ * 0) are the same quantities in the limit of fi — > 0. Dividing by (W)^^ 
in Eq. (|2.19|) insures that we remove the contribution to the Wilson loop not generated by the 



interaction with the shock wave |52| 



3. Static Solution 

As we discussed in the previous section, the calculation of the DIS cross-section at strong coupling 
under the AdS/CFT correspondence reduces to finding the classical trajectory of the string in the 
background of the shock- wave metric in Eq. fl2.13p . However, finding the time-dependent analytic 



solution to the Euler-Lagrange equations associated with the Nambu-Goto action in Eq. (|2.7|) in 
the background of the metric ( p,13[) is a very difficult mathematical problem. In this section we 
concentrate in a somewhat simplified problem, the static limit, which we expect to be a good 
approximation for DIS on a very large nucleus. 
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Following the standard procedures used in perturbative calculations || |8|, 25 1 we begin by 



relaxing the ultra-relativistic limit, allowing the nucleus to have a finite longitudinal extent in the 
x~ direction. This can be managed by "smearing" the delta function in Eq. ( |2.12| ): 



N 2 N 2 u 

(T__)= -^(ar)— > -JLao{x-)0(a-z-), (3.1) 



where 9 is the Heaviside step function and 



A A 1 ' 3 

2R-~ (3.2) 
p+ p+ 



is the width of the nucleus as seen from the dipole rest frame. (R is the radius of the nucleus.) After 
the replacement in Eq. Q3.1D the shock wave metric in Eq. ( |2.13| ) has to be modified accordingly, 
getting 



2 



ds 2 = \ 



-2dx + dx + —9(x )6(a — x ) z 4 dx 2 + dx\ + dz z 
a 



(3.3) 



Recalling that we considered the string to be located at x 3 = 0, Eqs. ( |3.1| ) and ( PT3| ) permit the 
following neat description of the temporal sequence of the collision. For t < there is no interaction 
between the string and the nucleus and the string lives in the empty AdSs space. At t — the 
string is hit by the front end of the nucleus. The interaction between the nucleus and the string 
continues until t = a \^2. At that time the nucleus has completely passed through the string, which 
returns to the empty AdS 5 space. We now concentrate on times t > 0. We will assume that the 
string reaches a stationary time-independent configuration in a short excitation time t e after being 
hit by the nucleus. Under that assumption, the string will remain in such stationary configuration 
during the remaining interaction time, i.e for t e < t < a\pl. Finally, at t > a \/2 the interaction 
ceases and the string is likely to return to its original vacuum configuration after a given relaxation 
time, t r . Very schematically, the contribution to the string action due to the interaction with the 
nucleus, 

S I NG = S NG (fi)-S NG (0), (3.4) 

can be written as 

S J NG = { j ' dt + J dt+ J dt) / dxC (3.5) 

where £ is the Lagrangian density. Under the additional assumption that both the transition and 
relaxation times are much smaller than the nuclear width, i.e. that 

a >> t e ~ t r , (3.6) 
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one can conclude that Sj^ G is dominated by the contribution of the static, time-independent solution 
(the second term in the r.h.s of Eq. ( |3.5|) ). Unfortunately, we were not able to find a reliable 
parametric estimate of neither t e nor t s in terms of the dimensionful scales of the problem, \x and a. 

Instead we notice that, as can be checked explicitly, the time-dependent equations of motion 
for a string right after the nucleus hit it (and right after the string leaves the nucleus) depend only 
on the ratio fi/a and on the boundary conditions at x — ±r/2, i.e., on r. (As we will see soon, the 
same applies to the static string equations of motion inside the shock wave, which do not "know" 
about the length of the shock wave.) As neither nja nor r depend on A (see Eqs. (|2.17|) and ( p. 2D ), 
we conclude that t e and t r , along with the Nambu-Goto action for the string, are independent of 



A. Hence the second integral on the right hand side of Eq. (3J3) is the only A-dependent term in 
Eq. ( p.5| ) and, as it grows with A, it is bound to dominate at large enough A. Therefore, in the 
limit of very large nucleus the static approximation will be justified. 

One could also justify the static approximation in a different way by slightly changing the 
physical problem. One could imagine that in DIS the virtual photon enters the shock wave, and 
then fluctuates into the qq pair. The qq pair travels through the shock wave, and, for the forward 
amplitude, it would quickly recombine back into a virtual photon before leaving the shock wave. 
Interaction of such a dipole with the shock wave has to be static, as the dipole would not "know" 
about the ends of the shock wave. Integrating over the possible splitting and recombination times 
would yield Eqs. (|2~20| ) and ( grp . 



In the static limit corresponding to scattering on a very large nucleus as described above all 
the relevant dynamical quantities become independent of time. We take the smeared large nucleus 
with the energy- momentum tensor described in Eq. fl3.1|) . The metric corresponding to such shock 
wave is given in Eq. (|3.3| ). To achieve the static limit one has to remove the product of theta- 
functions in Eq. (|3.3| ). Changing from light-cone coordinates to {t,x 3 ) coordinates and dropping 
the theta-functions in Eq. (|3.3|) we obtain 



ds l 



z 2 



(l - — z A ) dt 2 + f 1 + — z A ) (dx*) 2 -^zUtdx* + (dx l ) 2 + (dx 2 ) 2 + dz 2 ]. (3.7) 
\ 2a J \ 2a / a J 



Let us note that the shock wave metric in Eq. ( [3.7P , just like the metrics in Eqs. (|2.13j ) and (|3[ 
is an exact solution of the Einstein equations in empty AdS 5 space 

R^ + ^9^ = 0. (3.8) 

On the gauge theory side the metric in Eq. ( |3.7| ) corresponds to a uniform shock wave with an infinite 
longitudinal and transverse extent and with the only non-zero component of the energy-momentum 
tensor given by (T__) = ^ £ . 

Substituting the static metric Eq. (|3.7|) into the Nambu-Goto action ( |2.7|) and making use of 
the relation between the AdSs curvature radius L, the slope parameter a' and the string coupling 
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yields 



^ = VX (3.9) 



a 



a\f2 



rr/2 

S NG (fx) = I dt I dxC statlc , (3.10) 







r/2 



where 



£ *" = ^?V (1+i " ) ( 1 -5s* 4 ) < 3 - n) 

is the static Lagrangian density, z = z(x) is the string coordinate along the 5th dimension of 
the AdSs space (henceforth we shall drop the argument of z). The prime in Eq. ( |3.11D denotes a 
derivative with respect to x. 

We now have all the necessary ingredients to proceed to the calculation of the classical string 
configuration in the static limit. Before doing we note that 

£»ip*A». (3,2, 

The center of mass energy s of the collision is proportional to the product of p + and some scale 
characterizing the dipole. Indeed our dipole is made out of heavy quarks with large mass M, so it 
may appear natural to put s ~ p + M. However, later on we will renormalize the Nambu-Goto action 
that we will obtain by subtracting its contribution to the quark mass. Hence we will effectively 
remove M out of the problem. Hence the center of mass energy of the dipole would be equal to p + 
times some transverse momentum scale characterizing the dipole. We will put the latter scale to 
be proportional to the scale A of the nucleus, such that s ~ p + A. We thus "define" 

fa = *■ < 3 ' 13 > 
The Euler-Lagrange equation of motion associated with the Lagrangian Eq. ( EUTD reads 

d Q/^ s ^ a ^ c Q£static 

dx dz' dz ^ 

and straightforwardly leads to 

zz"(l-s 2 z 4 ) +2(z' 2 + l) = 0. (3.15) 
We need to solve Eq. ( 3.15Q with the boundary condition 

z(x = ±r/2) = (3.16) 



13 



which insures that the string ends at the boundary of the AdSs space. 

Using the first integral of Eq. ( |3.15| ) can be readily obtained. We get 



J2 



4 1 - s 2 z A 

1 _ s 2 z 4 

x 15 ^max 



- 1 



(3.17) 



where the constant of integration z max has the meaning of the maximum extent of the string in the z- 
direction, since z' = at z = z max (see Fig. |3|). z max is the fundamental parameter characterizing the 
static solution since, as we will show shortly, all the relevant physical quantities can be parametrized 
directly in terms of z max . 

As dx = dz/z', we write _ r /2 ft +r/2 



dx 



dz 



I / ( %max \ ^ 



l-s 2 2 4 



(3.18) 



- 1 



Picking the half of the string where z' > we integrate over 
x from — r/2 to x and, using Eq. (|3.16|) obtain 



x + 



dz 




I ( Zmax \ 



(3.19) 



Performing the integral in Eq. ( |3.19| ) yields 

-v3 



r 

x + - 



2 3 z 3 



s 2 Z A 

" max 



Fl 1 - 3 - 7 -^- 

2'4'4'zi.. 



Figure 3: A sketch of the string con- 
figuration solving Eq. ( |3.15| ) and given 



by Eq. (3.23) in the (x, z)-plane. 



(3.20) 



with F the hypergeometric function. The curve in Eq. ( |3.20| ) 
is sketched in Fig. [3|. 

From symmetry considerations z max must lie at the midpoint between the two string endpoints, 
i.e. at x — 0. Putting x = and z = z max in Eq. (|3.20|) yields the following relation between z max , 
the collision energy y/s and the dipole size r 



c r 



s 2 Z A 

" max i 



where we have defined a constant 



c = 



r 2 Q) 

(2tt) 3 / 2 ' 



(3.21) 



(3.22) 



The case considered by Maldacena in |62| for a static Wilson loop without a shock wave can be 
recovered from Eq. (|3.21| ) by putting s = in it, thus eliminating the shock wave. In that limit one 

gets Z max %max = ^0 
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Finally, with the help of Eq. ( |3.21|) , the classical solution for the open string in Eq. ( |3.20|) 

becomes 



r 

X + 2 



Q 7 3 



,13 7 z 4 

c o r F n> t; 7l 



2 7 4 ' 4 ' z 4 

max 



(3.23) 



for — (r/2) < x < with a mirror image obtained by replacing x —>■ —x in Eq. (|3~23|) forO < x < r/2. 
4. Evaluating the S- matrix 



We proceed by evaluating the Wilson loops in Eq. ( |2.19| ) using Eq. ( |2.10| ) with the classical string 
configurations found above. As we will see now, Eq. fl3.21| ) relating the string maximum to the 
collision energy and the dipole size has three different complex-valued solutions. As we shall see 
below, each of the solutions leads to a different behavior of the DIS cross-section. It turns out that 
the physically meaningful solution can be obtained in two ways which are described below. 

4.1 General Evaluation 



We first study the solutions of Eq. (|3.21|) for the string maximum, z max . Introducing the dimen- 
sionless parameters 



2 9 ' 
Cq r z 



4 4 2 

m = Cn r s , 



(4.1) 



Eq. ( |3.21D can be rewritten as the following cubic equation for £ 



mf d -£+ 1 = 0. 



(4.2) 



Obviously, for m 7^ Eq. ( f4.2|) has three different complex roots. Solving the cubic equation ( |4.2| ) 
we get 



1 



3mA 



+ A, 



(4.3) 



with 



A 



1 / 1 
+ 



1 



2 m V 4 m 2 21 m A 



1 1/3 


2irn 




exp 


1 

3 



(4.4) 



The index n in Eq. (|4.4j ) takes on values n = 0,1,2, corresponding to the three different Riemann 
sheets of the cubic root, and characterizes the three solutions of Eq. (|4.2| ). From here on we will 
stick to the convention that the cubic root is evaluated by taking the principal branch, while all 
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the information distinguishing one branch from another will be shown explicitly in terms of n. Also 
one writes, using Eq. (f4.1|), 



c r 



VI, 



(4.5) 



where we agree to take the principal branch of the square root. (One can show that the secondary 
branch does not lead to any new physical solutions.) The branch cut for all roots is taken along the 
negative real axis. 



Re [W(r=l GeV- 1 , s)] 



Im[z max (r=lGeV-\s)] 



n=0 
n=l 
n=2 



t S (GeV 2 ) 



a S (GeV 2 ) 



Figure 4: Real (left plot) and imaginary (right plot) parts of z max (r = lGeV _1 ,s) given by Eqs. ( |4.5[ ), 
( |4.4| ) and fl4.3|) plotted as a function of the collision energy s. The figure contains all three branches of the 
solution corresponding to n = (dashed line), n = 1 (dot-dashed line) and n = 2 (solid line). 




Figure 5: Real (left plot) and imaginary (right plot) parts of z max (r,s = lGeV 2 ) given by Eqs. ( |4.5[ ), 
flO]) and (£D|) plotted as a function of the dipole size r. Again the three branches of the solution are n = 



(dashed line), n = 1 (dot-dashed line) and n = 2 (solid line). 

The three different z max corresponding to the three solutions of Eq. ( f4.2|) are shown in Figs. ^ 
and [| as functions of energy s and the dipole transverse size r. As shown in Figs. [| and [|, and as 



- 16 - 



can be demonstrated by an analytical calculation, one of the three solutions of Eq. (O), the one 



corresponding to n = 1, is purely imaginary while the other two, corresponding to n = and n = 2, 
are real for small r and for small s, becoming complex and conjugate to each other for larger r and 
s once the condition of m > 4/27 is satisfied. 

Complex-valued solutions like these arise in the quasi-classical approximations to quantum 
mechanics (see e.g. ch. 131 in J76|). They correspond to motion in the classically forbidden region. 
Our physical interpretation of the complex-valued string coordinates is that the string is classically 
prohibited from going through the shock wave: such motion can only happen quasi-classically. 
Another way of thinking about this result is to observe that the saddle-point corresponding to the 
classical approximation does not have to lie in the domain of real-valued string coordinates only, 



and can in general be complex- valued. For similar cases in AdS / CFT framework see JF7 . 

In order to figure out which of the three branches leads to the right physical behavior of the 
dipole scattering amplitude we need to know the limiting behavior of z max asm->0 and m — > oo 
respectively. A straightforward expansion of Eqs. (f4.5p, ((Oj) and ( f4.3p yields 



Vi) , A i/ 3 I ex P Hf] > iorn = ( Vv^, forn = 

lim z max = ( ) < i, forn = 1 ; lim z max = < i/^/s, forn = 1 . (4.6) 

exp |i~J , forn = 2 ^ Co r, forn = 2 



From Eq. ( [4.6[ ) one can see right away that only the n = 2 branch of the solution maps onto the 
Maldacena result from [B2| in the limit of no shock wave, i.e. as m — » 0. One can also see this 



from Figs. § and [|. Thus if we were to require that the string coordinates map onto Maldacena's 
case for m — > 0, it would seem clear that one has to take the n = 2 branch of the solution. 
However, as we will see shortly, such branch does not give a unitary N(r, Y) at large r: in fact 
N(r, Y) becomes negative at large r for the n = 2 branch. Hence, to obtain a physically meaningful 
solution, we have to either abandon the requirement of mapping the string coordinate onto the 
Maldacena configuration H52| , or allow for the string to "jump" from one branch to another. The 
two possibilities are explored below. 

First let us calculate the Nambu-Goto action at the extremal solutions found above. Combining 



Eqs. ( |3.10| ) and (|3.11| ) and using dx = dz/z' we write 



SnoM = — / ~y/(l + ^)(l-^) (4.7) 

7T J Z Z 



where the factor of 2 comes from adding the two halves of the string. Using Eq. ( |3.17| ) in Eq. (|4.7| ) 
yields 



The integral in Eq. ( |4.8| ) has an ultraviolet divergence at z = due to infinite quark masses. 
To renormalize the result one has to subtract out the infinity p|J by removing contributions of 
individual quarks. The mass of a single quark in vacuum is calculated by attaching a straight-line 
string to the boundary of the AdS space, with the other end of the string going straight into the 



-r/2 



+r/2 



bulk. However, similar to the case of a heavy quark potential at finite temperature ||65|| , it is not 
clear here whether we should also subtract such an infinite straight string configuration, or let the 
string end somewhere. 

We would like to subtract out contributions of static 
straight open strings, "hanging" from the quark and the 
anti-quark into the AdS space stretching up to z = +00. 
However the shock wave metric in Eq. (|3.7|) appears to have 
an analogue of a horizon. Namely, for z — z^ with the 
"horizon" defined by 



1 



(4.9) 



V 



Figure 6: The string configuration 
which we subtract out to renormalize 
the Nambu-Goto action drawn in the 
(x, z)-plane. The dashed line indicates 
the horizon of the shock wave. 



the g a component of the metric becomes zero. Moreover, for 
z > Zh, no light-like trajectories along the z-axis (without 
moving in x 3 direction) are possible. Therefore the string 
world-sheet can not stretch into the z > z^ region. The 
situation is analogous to the finite temperature calculations 
of heavy quark potential, where one subtracts the contri- 
butions of strings stretching up to the black hole horizon 

pB| , 155] . Following |k| |55J we will subtract the straight-string configuration shown in Fig. [3] out 
of the action in Eq. ( f4.8|) . Our strings in Fig. |] stretch only up to the horizon z^. The action 
corresponding to two straight strings shown in Fig. ^| is obtained from Eq. ( f4.7|) by taking z 1 — > 00 



(4.10) 



limit in it and limiting the integration between < z < z^. Thus we will subtract 



— Zh 
2 A a f dz 



7T 



Vl - S 2 Z A 



from the action of Eq. 

Subtracting Eq. ( |4.10| ) out of the action in Eq. ( |4.8[ ) and integrating over z we obtain the 
renormalized action 



\f\ a \ Cq r 2 



71" C 



y/2 



7° 7. 



2 2 
— + — 

raax %h 



(4.11) 



The action ( |4.11| ), when substituted into Eqs. ( |2.20|) and ( p,19| ) would give us the scattering 
amplitude N(r,Y). 
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Knowing the action (4.11), and therefore the amplitude N(r,Y), we can now see which ones 



of the three above solutions give physically meaningful N(r,Y). Our criteria are simple: we want 
N(r, Y) — > when r — > (color transparency) and N(r, Y) — > 1 when r — > oo (the black disk 
limit). In addition N(r,Y) should be positive-definite for r > 0, N(r,Y) should be less than 1 at 
all r (unitarity) and N(r, Y) should be a monotonically increasing function of r for all r. 

A simple analysis using Eqs. ( [4.11] ), ( |2.19| ) and ( |2.19| ) rules out the n = 2 branch: as one can 
show, N(r, Y) given by that branch becomes negative at large r. n = branch appears to satisfy 
both the color transparency condition at small r and the black disk limit at large r. However, 
the amplitude N(r, Y) on that branch is not a monotonic function and, more importantly, is not 
unitary: in a certain range of energies/dipole sizes it oscillates and could become greater than 1 
violating unitarity. Hence this branch is also unphysical and should be discarded. This leaves us 
only with the n = 1 branch, which, as we will see shortly, gives a physically meaningful N(r, Y). 

The drawback of the n = 1 branch is that it does not map onto Maldacena's solution |)2J in 
the \i — > limit. Requiring that such mapping takes place leaves us with the n = 2 branch, which 
is unphysical at large r. However, there is another possibility, which we will outline below. One 
can take a superposition of two branches for z max : for small r (or s) one can use the n = 2 branch, 
while for larger r (or s) one can use the n = branch. As the n = and n = 2 branches of z max 
intersect at m — 4/27, as can be seen from Figs. f| and |5] and from Eq. (|4.4j ), it is possible that a 
transition from one branch to another happens at this point. As we will also show below, to get a 
physical result from the superposition of the two branches one has to slightly modify the rule for 
calculating the contribution of the Wilson lines to the S"-matrix ( [2.19 ). We will consider the single 
n = 1 branch solution first though. 

4.2 Single Branch 

As one can see from Eq. ( |4.6| ), and from a simple graphical analysis of Eq. ( |4.2| ), and also from 
Figs. |] and |^, z max given by the n = 1 branch of the solution is purely imaginary, with a positive 
imaginary part. As the coordinate x of the string is real, we infer from Eq. (|3.23|) that the complex 



phase of z should be the same as that of z max . Therefore, for n = 1 branch z would be purely 
imaginary. Denote the positive imaginary part of z max of n = 1 branch by p m ax such that 

Z max = i Pmaxi Pmax > 0. (4.12) 

Now, as z is purely imaginary, all ^-integrals, including that in Eq. ( |4. 1Q ) should run along the 
positive imaginary axis. Therefore for the n = 1 branch we should replace Zh — > % Zh in Eq. ( |4.11| ). 
We rewrite Eq. ( EOT ) as 



SJ33(Ai,n = l) 



\/A( 



7TC 



V2 



clr 2 2 2 



o 



+ 



Pvn&x PmcLX 



(4.13) 
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As one can explicitly check S t ^q{^ — > 0, n = 1) = 0. Using Eqs. ( |2.19| ) and (]2.20|) yields 

N(r, s) = 1 — exp 



y/Xt 



7TC 



v/2 



Car 2 



P 3 



+ 



2 



(4.14) 



where we replaced Y ~ In s with s in the argument of N. 

As one can explicitly check N(r, s) in Eq. (|4.14| ) satisfies all physical criteria for a meaningful 
solution: it goes to zero as r — > 0, it approaches 1 as r — > oo and it varies monotonically in between, 
always staying positive for r > and less than 1. The fact that the N(r,s) in Eq. ( [4. 14] ) is a 
monotonic function of r for all s is shown in detail in Appendix |A]. 

Therefore, Eq. ( |4.14 ) is our answer for the forward scattering amplitude of a dipole on a nucleus 
given by the n = 1 branch. We will refer to it as Solution A. To find p ma x one has to use Eqs. 
( |4.12| ), ( |4.5| ), ( f4.4|) , ( }4.3|) and ( f4.1|) for n = 1, which can be summarized as follows: 



with 



and m = cir 4 s 2 . 



Pn 
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1 



3mA 



n=l 



+ A n=1 



(4.15) 
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(4.16) 




s =10 <:/2 GeV 2 ,k=l,...,14 
Solution A 



Figure 7: Dipole scattering amplitude iV(r, s) from Eq. (|4.14| ) plotted as a function of the dipole size r for 
collision energies s = 10 fe / 2 GeV 2 , k = 1, . . . , 14 (from right to left) with A = 20, A 1 / 3 = 5 and A = 1 GeV. 
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N(r,s) from Eq. ( |4.14 ) is plotted in Fig. |^ as a function of r for a range of center of mass 
energies s. According to Eq. ( pT2|) we put 



A 1 ^ — 



(4.17) 



with A = 1 GeV. We also put A = 20 and A 1 / 3 = 5, which is not unrealistic for DIS on a gold 
nucleus. As one can see from Fig. [7], N increases with increasing energy, but the growth slows down 
as energy gets very high. 



& 2 (GeV ) 



Solution A 



10 4 10 r > 
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Figure 8: Saturation scale Q 2 as a function of the center-of-mass energy s corresponding to the same 
values of parameters as chosen for the solutions shown in Fig. 0. The saturation scale is obtained by 
imposing the condition N(r = 1/Q s ,s) = 0.5 on N(r,s) in Eq. ( [4.14| ). 

The saturation scale can be defined by requiring that N(r = 1/Q s ,s) = 0.5. It is plotted in 
Fig. | as a function of energy. The saturation scale also grows with energy, but at very high energy 
becomes a constant. It is interesting to note that fitting the region of linear growth of the logarithm 
of the saturation scale in Fig. |8] with a straight line one obtains 

O 2 s 
ln-^f « 0.9 + 0.47 In— (4.18) 
A 2 A 2 

However the growth stops approximately at s ~ 100 GeV 2 and Q 2 ~ 20 GeV 2 corresponding to 
XBj ~ Q 2 s /s ~ 0.2. For higher energies/smaller xsj the saturation scale is constant with energy/xsj- 



This behavior is qualitatively similar to the "saturation of saturation" conjectured in |78| using very 
different (small coupling) physical arguments. 

It would make more physical sense and facilitate the comparison to other approaches if we find 
the dependence of the forward amplitude N(r, s) and of the saturation scale Q s not only on s, but 
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also on xbj- The Bjorken-x variable in DIS is 



XBj 



Q 1 



s + Q 2 



(4.19) 



with Q the virtuality of the photon. Replacing Q ~ 1/r [79| we write 



1 



1 + sr 2 



such that for sr 2 C 1 



and for sr 2 ^> 1 



« 1 — s r 



(4.20) 



(4.21) 



(4.22) 



To better understand our results in Eq. (|4.14|) and in Figs. [7| and [8], let us study the asymptotics 
of N(r,s). At lower energies s and/or for smaller dipole sizes r we use the small-m asymptotics 
(see Eqs. (fO]) and (p~6|)) of the solution to write 



N(r, s) 



1 — exp 



r 2 s<l 
1 ™ 2 



V^Co 
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(4.23) 



Note that in this limit x B j ~ 1 — sr = o(l) and the conclusions derived from Eq. ( |4.23j ) will 
not apply to small-a;^ physics. The low-s and low-r asymptotics of the dipole is certainly energy- 
dependent, as can also be seen from Fig. [7]. One can use Eq. (|4.23| ) to extract the saturation scale: 
requiring that at r = 1/Q S the power of the exponent is approximately of order 1 yields 



Q 2 S ~ ^AA^A m < 1, 



(4.24) 



in agreement with the fit fl4.18|) of Fig. |] and the result of [|60[| . 

However, as m « s 2 r 4 , the condition of m 1 means that sr 2 C 1. We see that Eq. ( |4.23|) is 
valid only for large XBj- Rewriting Eq. ( |4.23| ) in terms of r and Xsj then yields 



N(r, s) 



r 2 s<l 



We see now that as a function of Xsj the saturation scale is 

Qs ~ v 7 ! - x B j A A 1 / 3 v^, x Bj < 1. 



(4.25) 



(4.26) 



The expression for the saturation scale in Eq. ( 4.26 ) is valid at large (order 1) x B j only. 
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At high energy and/or large dipole sizes corresponding to the small-x^j limit, Eq. ( 4. 14] ) in the 
m ^> 1 asymptotics gives 



N(r, s) 



l-exp<| - ^=rAA 1/3 \. (4.27) 
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Now x^j ~ l/(r 2 s) <C 1 (as sr 2 ^> 1) and we are in the sma\\-x Bj regime. The saturation scale 
can again be defined by requiring the power of the exponent to be of order 1. One then gets 

Q s ~ V\AA 1/3 . (4.28) 

The saturation scale is independent of energy s at high energies, in agreement with Fig. [S| It is also 
independent of xbj at small xsj- The saturation scale grows with the atomic number of the nucleus 
A as Q s ~ A 1 / 3 , which is a stronger growth than the perturbative QCD estimate of Q s ~ A 1 / 6 . 
The Q s ~ A 1 / 3 scaling was also observed in |60|| , though for the saturation scale more similar to 



the one given by our Eq. ( [4.24| ). The saturation scale similar to our Eq. ( |4.28| ) was found in |39] 



Finally we note that, while in both Eqs. ( |4.26|) and ( |4.28| ) the saturation scale grows proportional 
to \/\ at fixed xsj, the exact power of A depends on the energy- momentum of the shock wave used. 
As we have already mentioned above, the energy-momentum tensor of the gluon field at strong 
coupling has an extra power of v^A as compared to the energy-momentum tensor of the valence 
gluons C P- ID that we used above. Using the energy-momentum tensor of the gluon field in the shock 
wave metric would modify the powers of A in the above expressions for the saturation scales. In 
particular Eq. ( |4.28|) would then be modified to give Q s ~ A 3 / 4 . 



Let us stress that, while the solution presented in this Subsection and given by Eq. ( [4.14] ) is 
one of the two possible physical solutions for N(r, s), we believe that this solution is actually the 
correct one. As we will see below, the amplitude N(r, s) given by the n = 1 branch does not have 
discontinuous r-derivative, unlike N(r, s) given by the superposition of the other two branches. 
However, we do not have a solid physical argument to select the n = 1 solution over the other one. 

4.3 Superposition of Two Branches 



It is important to note that to derive our above result given by Eq. Q4.14 ) we had to abandon the 



condition that the string coordinates map smoothly back onto Maldacena's vacuum solution fl62 
when fi — > 0. Indeed the only branch of our solution satisfying such a condition, the n = 2 branch, 
would give a physically meaningless negative N(r, s) at large r. However, if one insists on smooth 
matching with the vacuum string configuration of jftllfl , one can construct a physically meaningful 
amplitude N(r, s) using a superposition of two branches: one can use the n = 2 branch of the 



solution given by Eqs. (|4.3|) , (|4. 1| ) and ( [4.4|) for m < 4/27 and the n = branch for m > 4/27. 



One can see from Eqs. ( f4.3| ), ( |4.1| ) and ( |4.4j ) that the values of z max given by those two branches 



are equal at m = 4/27, which insures smooth matching of the values of the Nambu-Goto action 
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between the two solutions. (The fact that the two branches intersect can also be seen in Figs. || 
and [5].) 

There is a subtlety here though. If one uses the Nambu-Goto action from Eq. ( |4.11| ) for the 
superposition of n = 2 and n = branches along with Eqs. (|2.19| ) and (|2.20| ) to find N(r, s), the 



resulting amplitude N(r, s) would still not be physical. In fact, in a range of values for s and r 
it still has oscillations, which make N(r, s) go above 1 and violate unitarity. These are the same 
oscillations that we have mentioned for the n = branch above. The origin of such oscillations is 
in the fact that z max for the n = branch has a non-zero real part (see Figs. (| and H), leading 
to a non-zero real part for the Nambu-Goto action, which, when substituted into Eq. Q2.19|) gives 
cosine-like oscillations. However, such oscillations can be removed by modifying the prescription 
for the calculation of the S"-matrix in Eq. (|2.19|) . One can replace Eq. Q2.19|) with 



S(r, Y) 



-Im[S NG M-S NG (»-+0)] 



(4.29) 



Indeed by explicitly keeping the imaginary part of the Nambu-Goto action in the exponent, as 
shown in Eq. ( |4.29| ), we would get rid of the real part and, therefore, of the oscillations described 
above. Also, the new prescription in Eq. fl4.29|) would not change our above conclusions in Sect. 
[4.2| as the action of the n = 1 branch is purely imaginary. However, we do not have a good physical 
justification of Eq. (^4.29|) in general. In the quasi-classical approximation to quantum mechanics 
one usually takes only the imaginary part of the action to estimate scattering amplitudes along the 



complex quasi-classical trajectories, as could seen from Eq. (131.14) in [[76]. Hence our prescription 
( (4.29| ) here is motivated by quantum mechanics and by the fact that it gives physically meaningful 
results here. It may be that for calculating the S-matrices, the prescription of p2~|| , which was 
originally proposed for heavy quark potentials, has to be modified to that of Eq. ( |4 . 2 9| ) . 
Substituting Eq. flPID into Eq. ( ^29|) , and using Eq. ( gg^ ) yields 



N(r, s) 



exp 



7TC 



V2 



Im 



2 



+ 



2 z r . 



(4.30) 



where we have modified the 2/zh term in the square brackets to reflect the fact that the complex 
phase of Zh (and of all other z's) should be the same as the complex phase of z max , as follows from 
Eq. ( |3.23|) and as we have noticed and used above in Sect. [12|. Zh is still given by Eq. (fl9|) above. 

Eq. ( |4.30| ) is the main result of this section. We will refer to it as Solution B. When using it one 
should take z max from the n = 2 branch of Eqs. ( |4.1| ), ( [4. 3| ) and ( |4.4| ) for m < 4/27 and from the 
n = branch of the same equations for m > 4/27. One can also notice that z max from both n = 2 
and n = branches is purely real for m < 4/27 (see Figs. |] and ^|). Hence in practice Eq. ( |4.30| ) 
can be used with the n = branch of z max for all m. This would eliminate the need for quantum 
corrections to justify the transition between the branches. However, the n = branch does not 
map onto Maldacena's vacuum solution p2[ : hence by keeping n = branch only we would also 
lose some of the justification for searching for solutions beyond the n = 1 branch described above. 
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Let us also note that, since z max and, therefore, the string coordinates for this solution are 
real for m < 4/27 and become complex for m > 4/27, one can interpret m = 4/27 as the point 
when, in a purely classical sense, the string would break due to high energy of the shock wave. This 



interpretation is similar to the finite temperature case Jj^, p% . 

The amplitude in Eq. Q4.30D is plotted in Fig. || for the same set of parameters as the amplitude 
in Fig. [7]. The striking feature of the amplitude in Fig. |9] is that N(r, s) = for a range of non-zero 



N(r) 




r(GeV~') 



Figure 9: Dipole scattering amplitude N(r, s) from Eq. ( 4.30| ) plotted as a function of the dipole size r for 
collision energies s = 10 fe / 2 GeV 2 , k = 1, • • • , 14 (from right to left) with A = 20, A 1 / 3 = 5 and A = 1 GeV. 



values of r. This means that non-zero transverse size qq dipoles below certain critical transverse size 
would not interact with the shock wave. While such behavior appears somewhat unphysical, one 
may interpret it by arguing that at large coupling most partons are located at small- x bj • As small 
r corresponds to large xsj (see Eq. (|4.20|) ), the small-size dipole does not resolve any partons, and 
thus does not interact. The exponential suppression of high-Q 2 (small-r) partons has been noticed 
before in ^3], [89[| . Such exponential suppression possibly translates into in our case. 

Eq. (|4.30|) can be used to find the saturation scale by requiring again that N{r = l/Q s , s) = 0.5. 

1 branch in Fig. |8|, the saturation 
Fitting the growing 



It is plotted in Fig. [Il] as a function of energy. Similar to the n 
scale also grows with energy, but at very high energy becomes a constant, 
part of the curve in Fig. [II] with a straight line yields 



In 



91 

A 2 



0.25 + 0.95 In 



A 2 ' 



(4.31) 



Again the growth of the saturation scale stops at a rather large values of XBj- 
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Qs (GeV 2 ) 



Solution B 



s (GeV ) 



Figure 10: Saturation scale Q 2 as a function of the center-of-mass energy s for the amplitude plotted in 
Fig. ^ and corresponding to the same values of parameters as chosen for the solutions shown in Fig. ||. The 
saturation scale is obtained by imposing the condition N(r = l/Q s ,s) = 0.5 on N(r,s) in Eq. (|4.30| ). 



Let us study the low and high energy asymptotics of Eq. ( |4.30| ). N(r,s) in Eq. ( |4.30| ) and in 
Fig. H becomes non-zero at m = 4/27 corresponding to 



4 



1/4 



27 J c y/s' 

Expanding N(r, s) from Eq. ([4 . 3 0|) for the n = branch in r — r* yields 



(4.32) 



N(r, s) 



1 — exp 



V / Aa2 1 /S 3 / 4 



r-r*<?;l/ v / s, 



r>r* 



IT C, 



1/2 



3V8 



(4.33) 



(We have not expanded the exponent here and in previous expansions as a may be large.) As 
r* ~ 1/a/s, we obtain ~ 1/r* 2 ~ s, as observed in Eq. ( (4.31| ). The Q 2 S ~ s scaling is somewhat 
unusual, as in terms of XBj it implies the existence of saturation Bjorken-x, some x B ^, such that 
the saturation sets in for XBj < x S Bj independent of Q 2 . However, as the Q 2 ~ s scaling is observed 
for large-^sj, we can not apply this conclusion for small-x^j physics, which we are studying using 
the Wilson loop. 

One can easily check using Eqs. ( Op and ( |4.30| ) that the s r 2 ^> 1 asymptotics of Eq. ( |4.30| ) is 



given by Eq. (|4.27|) . Thus the high energy asymptotics of both possible solutions discussed here is 
the same. The saturation scale at large s is also the same for n — 1 branch and for the superposition 
of n = 2 and n = branches. It is given by Eq. ([4.28 ) 
two possible solutions is at small r /low energy s. 



Therefore, the main difference between the 
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Let us finally note that, in DIS, for the Wilson loop to span the whole shock wave in the 
longitudinal direction, the interaction of the qq pair should be coherent over the longitudinal length 
of the whole shock wave. In other words, in the rest frame of the nucleus, the longitudinal coherence 
length of the qq pair should be larger than the diameter of the nucleus. The standard condition for 



this to happen in DIS is (see |Hl and references therein) 



1 



> 2R 



(4.34) 



2m N x B j 

with rriN the nucleon mass. Omitting factors of 2, putting ~ A and using xbj from Eq. ( |4.20| ) 
yields 



sr 2 > A 1/3 . 



(4.35) 



Violation of this condition does not imply a breakdown of our approximation. It would only mean 
that the Wilson loop does not stretch across the whole shock wave in the longitudinal direction, 
but only over a part of it. This would imply that, when the condition ( [4.35D is violated, A 1 / 3 in 
our above expressions for N(r, s) should be replaced by sr 2 . While this prescription would modify 
some of our l&rge-x Bj conclusions above, such as Eq. ( [4.23 ) along with the shape of the rise of the 
saturation scale in Fig. || our main small-x^ conclusion of constant Q s at high energy (|4.28|) would 
remain the same. 



5. Conclusions 

Let us restate our main results. We have studied DIS on a large nucleus in Af = 4 SYM using the 
AdS/CFT correspondence. We modeled the nucleus by an (infinitely) long shock wave. We argued 
that since for DIS in QCD the incoming virtual photon splits into a quark-anti-quark pair, which 
at high energies interacts with the nucleus eikonally (without recoil), the DIS cross section can be 



related to the expectation value of the Wilson loop, as shown in Eqs. (|2.1|) , fl2.2|) and (|2.3|) . We have 
calculated the expectation value of the Wilson loop by minimizing the area of the world-sheet for 
an open string connecting the quark and the anti-quark lines in the loop. We obtain two physically 
possible solutions resulting from different complex branches of string coordinates in the bulk. The 
two results for the forward scattering amplitude N(r,s) are given in Eqs. ( |4.14| ) and ( f4.30| ), which 



are the main results of the paper. They are plotted in Figs. [7] and |9|. The corresponding saturation 
scales are plotted as a function of energy in Figs. |8] and |10[ In both cases the main feature is that 
the saturation scale becomes independent of energy/Bjorken-x at sufficiently high energies/small 
values of xsj- 

We interpret this result as follows. At small gauge coupling, say in the BFKL evolution, 
each gluon emission is suppressed by a power of a s , but is enhanced by a power of Ins, which 
makes the resummations parameter a s Ins of the order of 1. High energy s is needed to generate 
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perturbative result 
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our large-coupling result 



In A 



QCD 



InQ 



Figure 11: A sketch of the saturation line in the (In Q 2 , In 1/xBj) plane. In the small coupling regime 
(Q 2 3> Aq CD ) the saturation scale grows according to CGC predictions. In the large coupling regime 
(Q 2 < Aq CD ) the saturation scale is a constant function of xsj in accordance with our results in this work. 



this enhancement. However, at large coupling the enhancement by logarithms (or by any other 
functions) of energy is not needed anymore, as A ^> 1. Hence one might conjecture that at large A 
partons fill up the whole phase space just due to sheer strength of the coupling. This way no energy 
enhancement is needed and increasing the energy should not change anything in this picture. This 
would lead to saturation scale being independent of energy. 

As our calculation was done for Af = 4 SYM, it is important to understand its implications 
for QCD. In an effort to do so we have plotted the saturation line in the (InQ 2 , Inl/xgj) plane 
in Fig. |ll|. The traditional physics of CGC happens in the small coupling regime of QCD for 
Q 2 ^> Aq CD , where Aqcd is the QCD confinement scale. The saturation scale there grows with l/xsj 
as some power Q 2 ~ (1/ ' XBj) 6as+ "' in the LLA, with ellipsis denoting the higher order corrections. 
The strong coupling limit of QCD is Q 2 < Aq CD . If one assumes that the large-coupling analysis 
carried out in this paper is applicable to the non-perturbative large- coupling Q 2 < Aq CD sector of 
QCD, then we can conclude that in this region Q 2 ~ (1/xbj) ■ This is pictured in Fig. 11. In 



the non-perturbative region we have simply sketched the curve from Figs. |8| and [10]: saturation 
scale starts out growing with decreasing XBj, but then levels off to a constant. As XBj decreases 
further, perturbative saturation mechanism would turn on and small transverse size partons would 
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be produced copiously due to lnl/xsj enhancement, leading again to the saturation scale growing 
as a positive power of 1/xgj. Indeed Fig. [11] is applicable only for a relatively small/dilute target, 
which would insure that there is a region of small xbj where Q 2 S is still in the non-perturbative 
region. For a very large nucleus or at very high energies perturbative QCD effects would dominate 
making Q 2 S perturbative. At the modern-day accelerators, such as HERA, Tevatron, RHIC, and 
LHC, one gets Q s ~ 1 4- 3 GeV for protons and nuclei [5B|. This scale may not be large enough 
to completely place us in the perturbative QCD region. Therefore it is likely that the physics at 
these machines lies in the transition region between the perturbative and non-perturbative regimes 



pictured in Fig. 11 



Addendum 

After the paper was already in print, we realized that by modifying certain mathematical conventions 
one could obtain the dipole scattering amplitude N(r) given in Fig. 9 from a single branch of the 
solution for z max as opposed to the superposition of two branches described in Sect. 4.3. A more 
conventional definition of the Nambu-Goto action in Minkowski space is different by a minus sign 
from our Eq. (2.7). Using only the n = 2 root for z max , denoted z 7 ^ 2 ., in such action along with Eq. 
(4.29) would lead to the dipole amplitude in Fig. 9. Alternatively, using the complex conjugate of 
the n = 2 root for z max , denoted z^^.*, in Eq. (4.30) would also give the dipole amplitude in Fig. 
9 coming from a single branch of the solution. This observation does not modify any of the results 
in the paper (and in particular in Sect. 4.3). 
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A. Monotonicity of N(r, s) as a function of r 

We wish to show that N(r, s) appearing in Fig. [7] is a monotonic function of r for all s, and not 
just for the values of s used in the plot. In particular, we claim that N(r, s) goes to zero when 
r —>■ and N(r, s) goes to 1 when r — > oo with a monotonic interpolation between the two extremes 
and without any local minima or maxima in-between. This would imply that N(r, s) is not only 
monotonic, but in addition is restricted between zero and one, a necessary condition arising from 
unitarity. 
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The minimum requirement we need to prove is that the derivative of the expression in the 
square bracket in the exponential of (|4.14j) is positive for any (positive) values of r and s. We 
should have in mind that p max = p ma x(r,s) [see Eqs. ( |4.12| ), (fO]-fO])]. We denote by E(r,s) the 
expression in the square brackets in the exponent of ( 4.1 4|) 



r 2 r 2 2 2 
E(r, s) ee 3 + - —-. (Al) 

PmaxVi^) Pmax 



We proceed by evaluating its derivative with respect to r 



dE _ f3c 2 r 2 + J2_\ %_ + 2c|r 



<9r \ o 4 p 2 7 <9r o 3 

\ rmax rmax / rrr 

At this point we recall that z maa; = ipmax- Eq. (|3.21|) then gives 



C r — S P max Pmax- (A3) 

Differentiating both sides of Eq. ( |A3| ) with respect to r we obtain 

dp r , 



c n r 



-fl ■ rmax y ~ rmax - ) q 

Solving Eq. dA4|) for and eliminating ^f^- in Eq. (|A2"D yields 

Pmaxd E Zpmax + Zclr 2 



Pmax (3s 2 pLx - 1) — ( A4 ) 



+ 2. (A5) 



c 2 r <9r 3s 2 p^ a:c -/4 
We now eliminate p^aa; i n the denominator of Eq. ( |A"5| ) using Eq. ( |A3|) . This yields 

Pmax dE _ 2p 2 max + 3c 2 r 2 



c 2 r dr 3 (c 2 r 2 + p 2 max ) - p 



max 



+ 2 = 1 > 0. (A6) 



This completes the proof of our statement. As can be easily shown (see also Fig. ||) p max > 0. 
Therefore we conclude that 

w>° < A7 » 



for all r > 0, which means that N(r, s) in Eq. ( 4.14 ) is a monotonically increasing function of r. 
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